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SECTION  I 


INTRODUCTION 

The  computer  code  for  the  stress  analysis  of  aircraft  tires  is 
designed  to  solve  the  following  problems: 

•  Inflation  of  a  lifted  but  unloaded  tire,  and 

§  Rotation  of  an  inflated  but  unloaded  tire. 

The  code  is  subdivided  into  six  overlays  within  the  framework  of 
dynamic  storage  allocation.  In  the  first  three  overlays,  the  input  data 
is  reduced  to  set  up  quantities  associated  with  the  geometrical  config¬ 
uration.  The  fourth  overlay  calculates  the  element  stiffness  and  load 
matrices,  which  are  assembled  In  the  fifth  overlay.  The  sixth  overlay 
contains  a  direct  equation  solver  with  one  right-hand  side. 

In  the  following  section,  the  structure  and  modulation  of  the  code 
will  be  described  in  detail. 

A.  Storage  Allocation  and  Input/Output 

Characteristics 

The  data  management  of  the  computer  code  incorporates  those  pri¬ 
mary  features  of  the  CDC  6600  system  which  are  necessary  for  the  effi¬ 
cient  flow  of  large  sets  of  information.  Information  storage  and  retriev¬ 
al  procedures  were  designed  to  minimize: 

•  Central  memory  required, 

§  Input/output  access  time,  and 

•  Program  maintenance  and  modifications. 

The  following  main  features  of  the  CDC  6600  system  were  used  to 
achieve  the  above  objectives: 
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•  Random  access  input/output  subroutines, 
e  Unblocked,  unbuffered  files,  and 
t  Blank  common. 

The  random  access  subroutines  are  library  input/output  routines, 
supported  by  CDC,  which  provide  the  capability  for  direct  storage  and  re¬ 
trieval  of  records  on  a  file,  as  opposed  to  sequential  files  which  re¬ 
quire  accessing  the  records  preceding  the  desired  one.  In  the  computer 
code,  these  mass  storage  routines  are  used  extensively  to  store: 

§  Input  data,  and 
a  Computed  data 

between  overlays,  which  allows  the  program  to 
a  Select, 
a  Input,  and 
0  Output 

only  those  data  which  are  necessary  in  the  particular  overlay  under  con¬ 
sideration. 

Unblocked,  unbuffered  files  are  used  to  store  the  intermediate 
data;  for  instance,  in  the  calculation  of  the  element  stiffness  and  ele¬ 
ment  load  matrix.  These  files  are  accessed  repeatedly  in  a  sequential 
manner.  These  unblocked,  unbuffered  files  are  efficient  for  reading  and 
writing  large  records,  since  the  information  is  read  directly  from  the 
disk  into  the  program  array  area.  Note  that  the  records  in  blocked  and 
buffered  files  would  first  be  read  into  an  intermediate  system  storage 
of  the  central  memory  area,  and  then  transferred  into  the  program  array 


area. 


The  proper  use  of  blank  common  allows  the  code  to  have  a  general 
work  area  available,  whose  length  depends  only  on  the  field  length  declared 
on  the  job  card.  This  area  is  dynamically  divided  among  the  arrays  needed 
in  executing  the  current  overlay. 

Thus,  the  central  memory  disposition  of  the  codes  has  the  following 
structure: 

word  Iv _ 

®  system  area 

labelled  common 

program  instruction,  variables, 
arrays,  and  file  work  area 

blank  coimion 

field  length) 

The  above  construction  allows  the  user  to  specify  a  field  length  tailored 
to  the  data  size. 

In  the  current  code,  each  overlay  determines  the  length  of  the  ar¬ 
rays  used  and  stacks  them  nose-to-tail  in  blank  common. 

B.  Fortran  Extended  Code 

The  code  is  originally  developed  under  Scope  3.3  for  the  FTN  com¬ 
piler.  It  contains  both  Fortran  and  Compass  routines. 

C.  Library  Routines 

Besides  the  standardized  library  routines,  the  codes  employ  the 
following  special  features  of  the  CDC-6600  library: 


St  word  (  =  declared 
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•  BUFFER  IN, 

•  BUFFER  OUT, 
e  READMS,  and 

•  WRITMS. 

D.  Assembly  Language  Subroutines 

The  CDC-6600  assembly  language.  Compass  (COMPrehensi ve  Assembly 
language), is  particularly  suited  to  substantially  reducing  the  computa¬ 
tion  time  of  looped  operations.  The  improvements  are  realized  by: 
e  More  efficient  retrieval  of  array  elements, 

•  Overlapping  of  data  storage  and  retrieval  from  central  memory 
with  multiplication,  addition,  and  subtraction,  and 

9  Efficient  use  of  the  instruction  stack  which  holds  seven 
words  (up  to  28  instructions)  in  the  central  processor. 

Well -coded  Compass  routines  will  execute  computational  do-loops 

from  five  to  six  times  faster  than  normal  Fortran  IV  on  the  CDC-6600. 

In  the  present  code,  six  of  the  heavily  used  matrix  manipulation 

subroutines  are  written  in  Compass: 

e  MATMPY, 

•  MATADD, 

•  MATSMP, 

9  INPRDS, 

0  VECMAT,  and 
9  EMULT. 

There  are  six  special  Compass  subroutines  to  perform  tasks  which 
standard  Fortran  IV  is  not  designed  to  handle.  These  are  listed  below: 


ifil'i  i  fifilliiiillifi'iiiiifiillhili  Mi'tli~'iii  t 


•  KFL. 

•  SSZERO, 

t  MSTG, 

0  GET, 

•  PUT,  and 

•  STRMOV. 

KFL  is  a  Compass  subroutine  which  retrieves  the  field  length  request- 
ed  by  the  job.  This  Infonsatlon  allows  the  program  to  use  all  central 
memory  available.  Furthermore,  for  each  data  set,  a  minimum  field  length 
requirement  may  be  tailored. 

SSZERO  is  designed  to  set  array  values  to  zero  during  execution  in 
a  minimum  amount  of  time.  The  routine  Is  used  extensively  throughout  the 

code. 

MSTG,  GET,  PUT,  and  STRMOV  are  Compass  subroutines  which  perform 
character  and  string  manipulation.  They  are  used  in  dynamic  storage  allo¬ 
cation  and  by  the  free  field  input  reader  which  reduces  considerably 
the  time  required  by  the  user  to  enter  and  debug  his  input  data. 


SECTION  II 

PROGRAM  ORGANIZATION 


The  program  is  organized  according  to  the  following  problem 

types : 

•  Inflation,  and 
e  Inflation  and  rotation. 

The  general  data  flow  is  shown  below,  where  the  numbers  refer 
to  the  overlays  under  consideration. 
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SECTION  III 


OVERLAY  DESCRIPTION 

A.  Overlay  (KTIRE,  0.0) 

This  overlay  controls  the  general  data  flow  as  described  in  Section 
II.  It  performs  the  initialization  of  labeled  common  blocks,  opens  random 
access  mass  storage  files,  and  facilitates  the  storage  requirements  within 
the  framework  of  dynamic  storage  allocation.  This  overlay  also  contains 
various  utility  programs  and  assembly  language  routines  for  vector  and  char¬ 
acter  manipulations. 

B.  Overlay  (KTIRE,  1,0) 

All  the  input  data  are  read  from  cards  in  this  overlay  and  then 
they  are  placed  on  random  access  mass  storage  files.  The  input  data  are 
checked  for  logical  errors  which  are  summarized  at  the  end  of  the  data 
processing  phase,  using  the  subprograms  RANGE,  WRDCHK,  and  COMPCHK.  The 
primary  control  parameters  are  also  set  up  here  in  the  labelled  common 
blocks  SIZE  and  CONTACT. 

C.  Overlay  (KTIRE,  2,0) 

This  overlay  performs  the  preliminary  nodal  calculations  such  as 
t  Cartesian  and  curvel inear  coordinates, 

•  Surface  vectors  of  the  undeformed  reference  surface,  and 

•  Cord  angle  distribution  along  the  meridian. 

The  resulting  data  are  placed  on  random  access  mass  storage  files. 
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D.  Overlay  (KTIRE,  3,0) 

The  intrinsic  reference  element  properties,  such  as 
e  Element  area, 

0  Element  centroid, 
c  Local  unit  vectors, 

•  Local  element  vertex  coordinates,  and 
e  Average  cord  angles, 

are  generated  here  and  then  placed  on  random  access  mass  storage  files. 

E.  Overlay  (KTIRE,  4,0) 

This  overlay  sets  up  the  element  stiffness  matrix  and  load  vector 
according  to  the  finite  element  formulation  outlined  in  [1], 

F.  Overlay  (KTIRE,  5,0) 

This  overlay  constructs  the  structure  stiffness  matrix  for  all  the 
problems  under  consideration,  and  the  structure  load  vectors  for  infla¬ 
tion  and  rotation  problems. 

G.  Overlay  (KTIRE,  6,0) 

A  standard  gaussian  elimination  routine  is  contained  in  this  over¬ 
lay. 
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SECTION  IV 


FILES  AND  COMMON  BLOCKS 

The  program  control  variables  are  transmitted  via  labelled  com¬ 
mon  blocks,  while  the  relevant  input  or  calculated  data  are  placed  on 
sequential  or  random  access  mass  storage  files. 

A.  Common  Blocks 

BCINDEX  contains  the  information  for  dynamic  storage  allocation: 

•  NOPOS  or  number  of  positions  available  for  allocation, 

•  KSPACE  or  space  available  for  allocation,  and 

•  INDEX  or  pointers  to  currently  defined  arrays. 

CONTACT  consists  of  the  information  for  the  iterative  contact  alogorithm: 

•  NORING  or  number  of  contact  rings, 

•  KRING  or  current  ring  number, 

•  NODMAX  or  maximum  number  of  nodes  in  a  ring,  and 

•  KSZAX  or  total  number  of  candidate  contact  nodes. 

ERROR  accounts  for  all  logical  input  data  errors  and  includes: 
e  NERR  or  running  count  of  input  errors, 
e  NERRS  or  record  numbers  of  those  with  errors, 

§  NERLIM  or  maximum  number  of  errors  to  be  counted,  and 

•  KERR  or  indicator  as  to  whether  current  block  has  errors. 

FILES  contains  the  names  of  all  data  files  defined  by  Fortran  IV 
Hollerith  form: 
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INDIA 


consists  of  information  about  the  last  record  read  by  the  input 


reader: 

•  NWRD  or  number  of  items  present, 

t  ITYP  or  type  of  each  item  present, 

•  NREC  or  record  number, 

§  NCRD  or  card  number,  and 

•  DTA  or  value  of  each  item  in  the  record. 

MATSIZ  consists  of  the  structure  stiffness  matrix  characteristics: 

•  NIIMBK  or  number  of  blocks, 

•  NBKSI  or  block  size, 

•  NMIQ  or  bandwidth  including  right-hand  side, 

0  NPB  or  nodal  points  per  blocks, 

•  NEQ  or  number  of  equations  per  blocks, 
t  NMAX  or  total  number  of  equations, 

«  NORHS  or  number  of  right-hand  sides,  and 

•  NRBKSI  or  the  size  of  the  flexibility  matrix  block. 

PRINTS  controls  all  optional  intermediate  printing: 

•  KPRINT  or  an  array  indicating  which  intermediate  values  the 
user  wants  printed,  and 

e  LINLIM  or  the  maximum  number  of  single  spaced  lines  per  page. 
RECORD  contains  the  variable  names  for  all  named  random  access  records 
RETRIV  consists  of  information  to  determine  the  blocksize  of  the  struc 
ture  stiffness  matrix. 
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SIZE  is  comprised  of  the  input  control  parameters: 

•  NOEL  or  number  of  elements, 

•  NNODE  or  number  of  nodes, 

f  RADIUS  or  radius  of  meridian  reference  curve  rotation, 

§  NRHO  or  number  of  fitting  coefficients  for  the  meridian  ref¬ 
erence  curve, 

•  NPRHO  or  number  of  data  points  for  the  meridian  reference  curve, 

•  MAXLAY  or  maximum  number  of  layers, 

•  GREEN  or  green  angle, 

•  SPEED  or  rotational  speed,  and 

•  INCR  or  number  of  increments  for  nonlinear  inflation. 

B.  Random  Access  Files 

Following  in  Table  1  is  a  description  of  all  records  of  all  random 
access  files.  The  fourth  table  entry  indicates  whether  the  information 
is  input  or  calculated  by  the  code. 
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Table  1. 


Random  Access  Files 


FILE  NAME 

RECORD 

NAME 

CONTENTS 

I/C 

SIZE 

KCTMAT 

KCTC 

contact  data 

I 

NORING  * 
(NODMAX  +2) 

KCTMAT 

KCTC2 

x-coordinate 

array-contact 

case 

C 

KSZAX 

KCTMAT 

KCTC  3 

contact  case 

ALPHA  array 

C 

maximum 

KSZAX 

KCTMAT 

KCTC5 

total  nodal 
load  vector- 
contact  case 

C 

maximum 

KSZAX 

KHMAT 

numbered 
one/ layer 

decomposition  of 
the  HB  matrix 

C 

maximum 

378 

KLADAT 

numbered 

one/element 

element  layer 
data 

I 

maximum 
MAXLAY  *10 

KONDAT 

KAR 

elements  areas 

C 

NOEL  *1 

KONDAT 

KAV 

average  cord 
angles 

c 

NOEL  *1 

KONDAT 

KBETA 

curvefit  betas 

I 

NPBETA  *2 

KONDAT 

KCA 

cord  angles 

c 

NNODE  *1 

KONDAT 

KCAR 

cartesian 

coordinates 

c 

NNODE  *3 

KONDAT 

KCARC 

cartesian  data 

I 

NNODE  *6 

KONDAT 

KCEN 

centroi d 

c 

NOEL  *3 

KONDAT 

KCLC 

curvel inear 
coordinates 

I 

NNODE  *2 

KONDAT 

KDIS 

displacements 

I 

NNODE  *6 
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Table  1  (con't) 


FILE  NAME 

—imm - 

NAME 

CONTENTS 

I/C 

SIZE 

KONDAT 

KELN 

elements'  nodes 

I 

NOEL  *3 

KONDAT 

KFOR 

forces 

I 

NNODE  *3 

KONDAT 

KINCR 

nonlinear  incre¬ 
ment  data 

I 

INCR  *2 

KONDAT 

KIPV 

in-plane  vertex 
coordinates 

C 

NOEL  *6 

KONDAT 

KLUV 

local  unit  vectors 

C 

NOEL  *9 

KONDAT 

KRHO 

curvefit  Rhos 

I 

NPRHO  *2 

KONDAT 

KSV 

surface  vectors 

C 

NNODE  *9 

KPMAR 

numbered 

one/element 

upper  half  element 
load  and  stiff¬ 
ness  matrix 

C 

135  words 

KRMAT 

numbered 

one/layer 

H"*  matrix  after 
forward  substi¬ 
tution 

C 

maximum 
(NOLAY  *21-7) 
*21 

KSLMAT 

KSLV 

solution  vector 

C 

NNODE  *5 

KSTFIL 

numbered 

one/merged 

block 

output  merged  K 
matrices 

C 

maximum 

NEQ  * 

(KBAND 
+  KSZAX) 

2 

numbered 

one/block 

SOLVMOR  solution 

C 

maximum 

NEQ  * KSZAX 

SECTION  V 


INPUT  DATA  PREPARATION 

There  are  no  inherent  units  assigned  to  any  of  the  variables.  Thus, 
the  user  is  free  to  use  any  desired  set;  however,  once  defined,  the  units 

must  be  consistent  for  all  parameters. 

Format-free  processing  is  employed  for  both  floating  point  and 

integer  variables. 

Numeric  specifications  must  be  separated  from  each  other  by  at 
least  one  blank  column. 

'  A.  Title  Cards 

No  restriction  is  imposed  on  the  number  of  title  cards.  Each  of 
these  cards  must  start  with  a  slash  in  column  1.  Example: 

/  DATA  SET  NUMBER  1 
/  PREPARED  10/1/72 
/  AIRCRAFT  TIRE 
/ 

B.  Control  Cards 

A  name  is  assigned  to  each  of  the  control  variables.  The  name  is 
followed  by  the  appropriate  numerical  specification.  The  last  character 
on  each  of  these  cards  must  be  the  slash  symbol,  which  is  separated  from 
the  words  or  numeric  specifications  by  any  number  of  blank  columns. 

These  control  cards  follow  the  schedule  below: 

BEGIN  DATA  INPUT  / 

BEGIN  CONTROL  PARAMETERS  / 
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NODES  N  / 

ELEMENTS  E  / 

RADIUS  R  / 

NRHO  M  / 

NPRHO  N  / 

LAYERS  L  / 

GREEN  a  / 

SPEED  w  / 

INCREMENTS  IN  / 

END  CONTROL  PARAMETERS  / 

/ 

/ 

The  above  numeric  control  variables  are  defined  as  follows: 

N  =  Number  of  nodes  (integer), 

E  =  Number  of  elements  (integer), 

R  =  Rotation  radiu^ 

M  =  Number  of  curvefitting  parameters  for  the  median  section  (integer), 

N  =  Number  of  prescribed  data  points  for  the  meridian  section  (integer), 
L  =  Maximum  layer  number  (integer), 
a  =  Green  angle, 
u)  =  Rotational  speed ,  and 

IN  =  Number  of  increments  for  nonlinear  inflation. 

C.  Nodal  Data 

The  node  number  and  the  corresponding  curvel inear  coordinates  are 
specified  here: 
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BEGIN  NODAL  DATA  / 

i  e  <i) 

• 

END  NODAL  DATA  / 

/ 

/ 

The  puineric  veriables  i  i  snd  4’ 
i  =  Node  number  (integer), 

0=  Parallel  in  radians  (floating  point),  and 
(j)  =  Meridian  in  radians  (floating  point). 

D.  Element  Data 

There  are  three  nodes  associated  with  each  surface  element.  The 
node  assignments  must  follow  the  right-hand  rule  according  to  the  outward 
normal  direction.  Furthermore,  the  first  node  number  must  be  the  smallest 
Within  this  record,  the  elastic  constants  are  also  specified  for 
each  layer  within  the  element  in  a  principal  reference  frame.  Thus, 

BEGIN  ELEMENT  DATA  / 

E  j  N^  N2  N3  N^  / 

L  1  t  b  Sr  f  / 


L  k  t  r  Ep  Vr  E^  b  Sr  f  / 

END  ELEMENT  DATA  / 
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The  numeric  characters  in  the  above  records  are: 
j  =  Element  number  (integer), 

=  Node  numbers,  assigned  according  to  the  right-hand  rule  dic¬ 
tated  by  the  outward  normal  (integer)  (N-j  must  be  the  smallest), 

=  Total  number  of  layers  for  this  element, 

k  =  Layer  number  (integer), 

t  =  Layer  thickness, 

r  ■  Cord  versus  matrix  area  fraction  per  inch, 

E|^  =  Matrix  Young's  modulus, 

Vp  =  Matrix  Poisson  ratio, 

E^  =  Cord  Young's  modulus, 

=  Cord  Poisson  ratio, 
b  =  bias  multiplier, 

Sp  =  Rubber  mass  density, 

=  Cord  mass  density,  and 
f  =  Dummy  variable  for  future  use  (f  =  0). 

For  identification  purposes,  the  letter  E  must  precede  the  nodal 
information  and  the  letter  L  must  precede  the  layer  flexibility  informa¬ 
tion. 

E.  Prescribed  Forces 

At  a  given  node,  one  may  specify  three  external  force  components. 

The  input  records  read  as 
BEGIN  FORCE  DATA 
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END  FORCE  DATA 

/ 

/ 

where 

k  =  Node  number  (integer),  and 
=  Prescribed  force  components  (floating  point). 

F.  Prescribed  Displacements 

At  each  node,  three  rectilinear  and  two  rotational  displacement  com 
ponents  may  be  restrained  (i.e.,  equal  to  zero).  These  components  are 
referred  to  as  the  base  vectors  of  the  undeformed  reference  and  labelled 
according  to  the  following  schedule  for  a  rotational ly  symmetric  problem: 

=  Component  along  the  parallel, 
qg  =  Component  along  the  meridian, 
q^  =  Component  along  the  normal, 
q^  =  Rotation  component  along  the  meridian,  and 
qg  =  Rotation  component  along  the  parallel. 

The  input  records  are  exhibited  as  follows: 

BEGIN  DISPLACEMENT  DATA  / 
k  m  N"!  N2  . . .  N|y|  / 

END  DISPLACEMENT  DATA  / 

where 

k  =  Node  number  (integer), 

M  =  Total  number  of  displacements  restrained  at  the  node  under 
consideration  (integer),  and 


MpM2..  =  Number  of  displacements  being  restrained  (integer). 


Curvefitting  the  Meridian  Reference  Surface 


The  cartesian  coordinates  of  the  meridian  section  are  defined  in 
this  record.  Thus, 

BEGIN  CURVEFIT  RHOS  / 
i  X2  / 

END  CURVEFIT  RHOS  / 

/ 

/ 

where 

i  =  Sequence  number  of  data  points  (i  =  1 ,2,. . .NPRHO) (integer),  and 
Xi,X2  =  Cartesian  coordinates  (floating  point). 

H.  Increment  Data 

In  this  record,  the  increment  numbers  with  the  corresponding  incre¬ 
mental  pressure  are  defined  as  follows: 

BEGIN  INCREMENT  DATA  / 

I  p,  / 
i  p,  / 

IN  P,N  / 

END  INCREMENT  DATA  / 

where 

i  =  Increment  number,  and 
P^.  =  Incremental  pressure. 

I.  Print  Options 

The  output  information  is  governed  by  the  following  control  cards: 
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BEGIN  PRINT  OPTIONS  / 

ALL  / 

CONTROL  PARAMETERS  / 

NODAL  DATA  / 

ELEMENT  DATA  / 

CURVEFIT  RHOS  / 

LOCAL  UNIT  VECTORS  / 

NODAL  OUTPUT  TABLE  / 

INCREMENT  DATA  / 

CONTACT  DATA  / 

END  PRINT  OPTIONS  / 

/ 

/ 

END  DATA  INPUT  / 

J.  Commenting  the  Input  Records 

Following  the  slash  on  the  input  data  cards,  comments  may  be  inserted. 
These  comnents  may  be  continued  on  any  number  of  cards,  having  a  slash  for 
the  first  character.  Thus, 


BEGIN  NODAL  DATA  / 

1  0.5  0.8  -0.6  /  UPDATED,  NOV.  11,  1971 

2  0.6  1.0  -0.5  /  COWENTS 

/  COMMENTS 


/  COMMENTS 
/  COMMENTS 
3  0.8  2.0  0.5  / 


BEGIN  ELEMENT  DATA  / 

1  2  3  10  /  C0^^1ENTS 

/  COMMENTS 

/  COWENTS 
/  COMMENTS 

2  4  5  11  / 


K.  Input  Table  Summary 

The  input  data  card  set-up  is  suimnarized  in  this  section. 
/  DATA  SET  NUMBER  1 
/  PREPARED  10/1/72 
/  AIRCRAFT  TIRE 
/ 

/ 

BEGIN  DATA  INPUT  / 

BEGIN  CONTROL  PARAMETERS  / 

NODES  N  / 
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ELEMENTS  E  / 

RADIUS  R  / 

NRHO  M  / 

NPRHO  K  / 

LAYERS  L  / 

END  CONTROL  PARAMETERS  / 

/ 

/ 

BEGIN  NODAL  DATA 
i  0  (fi  / 

END  NODAL  DATA  / 

/ 

/ 

BEGIN  ELEMENT  DATA  / 

E  j  N^  N2  N3  N4  / 

L  k  t  r  Er  Vr  Ec  Vc  b  S 
END  ELEMENT  DATA  / 

/ 

/ 

BEGIN  FORCE  DATA 

k  Pi  P2  P3  / 

END  FORCE  DATA  / 

/ 

/ 

BEGIN  DISPLACEMENT  DATA  / 
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k  M  N2  N3  / 

END  DISPLACEMENT  DATA  / 

/ 

/ 

BEGIN  CURVEFIT  RHOS  / 

i  Xg  / 

END  CURVEFIT  RHOS  / 

/ 

/ 

BEGIN  CURVEFIT  BETAS  / 
1x^X3/ 

END  CURVEFIT  BETAS  / 

/ 

/ 

BEGIN  INCREMENT  DATA  / 
I  P  / 

END  INCREMENT  DATA  / 

/ 

/ 

BEGIN  PRINT  OPTIONS  / 
ALL  / 

CONTROL  PARAMETERS  / 
NODAL  DATA  / 

ELEMENT  DATA  / 

CURVEFIT  RHOS  / 
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LOCAL  UNIT  VECTORS  / 

NODAL  OUTPUT  TABLE  / 

ELEMENT  OUTPUT  TABLE  / 

END  PRINT  OPTIONS  / 

/ 

/ 

END  DATA  INPUT 

For  future  extensions,  the  code  angle  variation  may  be  described  by  exper¬ 
imental  data  points.  This  phase  of  the  code  is  not  yet  implemented,  how¬ 
ever,  as  the  corresponding  input  data  must  be  present.  Thus,  augment  the 
control  parameter  block  by 
NBETA  2  / 

NPBETA  2  / 
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SECTION  VI 
OUTPUT  DESCRIPTION 

A.  Output  Data 

During  the  data  processing  phase,  associated  with  geometrical  charac 
terization,  the  user  may  exercise  the  print  options  described  in  Section  V. 
The  corresponding  output  information  contains  the  following  records: 

•  Input  data  -  the  user's  input  is  listed  -  to  each  card  image  a 
record  number  and  a  card  sequence  number  are  assigned  for  error 
detection  purposes, 

•  Control  parameters, 

•  Curvel inear  coordinates  of  the  nodes, 

•  Element  and  layer  data, 

•  Cartesian  coordinates  of  the  data  points  of  the  reference 
meridian, 

•  Increment  data  for  nonlinear  inflation, 

•  Modal  output  table, 

•  Element  output  table,  and 

•  Local  unit  vectors. 

For  each  incremental  step,  the  generalized  nodal  displacements  are 
tabulated  as  follows; 

SOLUTION  VECTOR  DUE 
Node  U1  U2 

1  *  * 

2  *  * 
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TO  INFLATION 

U3  SI  S2 

*  *  * 

4r  ★  ★ 


where 


U1  =  displacement  in  the  parallel  direction, 

U2  =  displacement  in  the  meridian  direction, 

U3  =  displacement  along  the  outward  normal  to  the  reference  surface, 

51  =  rotation  with  respect  to  the  base  vector  along  the  meridian 

surface,  and 

52  =  rotation  with  respect  to  the  base  vector  along  the  parallel 

direction. 

The  above  information  is  then  followed  by  the  detailed  description  of  the 
state  of  stress  associated  with  each  element. 

ELEMENT  NUMBER 
n 

LOCAL  DISPLACEMENTS 
U  V  W  RX  RY 

*  *  ★  *  * 

*  *  *  *  ★ 

*  *  *  *  * 

where 

U  =  displacement  vector  in  the  gi  direction, 

V  =  displacement  vector  in  the  ^2  direction, 

W  =  displacement  vector  in  the  direction, 

RX  =  rotation  with  respect  to  g2,  and 
RY  =  rotation  with  respect  to  g’-j, 

where  ^1 ,  g2,  and  ^re  local  unit  vectors  described  in  Volume  I. 
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This  information  is  then  followed  by  the  list  of  stress  resultants 
and  reference  surface  strains  and  the  stress  distribution  across  the  thick¬ 
ness  as  follows; 


STRESS  RESULTANTS 

Nil  N22  N12  Mil  M22  M12  Q1  Q2 

******** 


REFERENCE  SURFACE  STRAINS 

Ell  E22  El  2  Kll  K22  K12  El  3  E23 

*****  *** 


STRESS  DISTRIBUTION 


Layer  Number 

Sll 

S22 

SI  2 

SI  3 

S23 

1 

★ 

* 

* 

* 

* 

2 

★ 

* 

* 

★ 

* 

where  the  stress  resultants  and  reference  surface  strains  are  calculated 
at  the  centroid  of  the  reference  surface  element.  The  corresponding  layer 
stresses  are  calculated  at  the  midpoints  of  the  layers  in  the  thickness 
direction. 

The  same  type  of  information  is  then  repeated  for  the  rotation  if 
requested  by  the  SPEED  control  card. 
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B.  Error  Exits 

Extensive  input  error  checks  are  provided  in  the  data  preparation 
phase.  Each  input  card  is  traced  according  to  its  internal  sequence  num¬ 
ber  in  the  input  deck.  For  cross  reference,  appropriate  record  numbers 
are  also  internally  assigned  to  the  input  cards.  Some  of  the  most  frequent 
error  types  are  listed  below: 


RECORD  1 )  BEGIN  DATE  INPUT  /  MISSPELLED  DATA 
—ERROR  10010  INPUTD  -  FIRST  DATA  CARD  MUST  CONTAIN  BEGIN 
DATA  INPUT  -  PROGRAM  TERMINATES 

RECORD  2)  BEGIN  DATA  INPUT  (FORGOT  SLASH) 

—ERROR  10020  INPUTD  -  BLOCK  NAME  MUST  CONTAIN  EXACTLY  THREE 
WORDS  -  NEXT  INTENDED  RECORD  IS  IGNORED  AND  PROGRAM  TERMINATES 

RECORD  3)  NODES  8.  /  DECIMAL  POINT  ..CARD  7 
—  ERROR  10040  INPUTD  -  THE  INPUT  VALUE  FOR  ELEMENTS  MUST  BE 
TYPE  INTEGER 

RECORD  4)  ELEMENT  12.  /  DECIMAL  POINT  ..CARD  8 
—ERROR  10040  INPUTD  -  THE  INPUT  VALUE  FOR  ELEMENTS  MUST  BE 
TYPE  INTEGER 

ERROR  SUMMARY  -  THERE  ARE  (IS)  TWO  ERROR  IN  INPUT  DATA  -  THE 
RECORD  NUMBERS  AT  WHICH  ERRORS  ARE  DETECTED  ARE  PRINTED 
BELOW  -  ONLY  THE  FIRST  100  DETECTED  ERRORS  ARE  LISTED 

3  4 


Similar  logical  errors  are  detected  and  summarized  at  the  end  of 
the  program  so  that  the  user  may  correct  several  input  errors  in  the  same 
run. 


C.  Timing  and  Field  Length 

The  field  length  requirements  may  easily  be  established  based  on 
the  available  blank  common  in  the  overlays,  whose  lengths  are  listed  below 
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(Scope  3.3): 


OVERLAY  (KTIRE,0,0) 

100  -  25116  (OCTAL) 

OVERLAY  (KTIRE,1,0) 

25116  -  37220 

OVERLAY  (KTIRE,2,0) 

25116  -  30400 

OVERLAY  (KT IRE, 3,0) 

25116  -  27257 

OVERLAY  (KTIRE,4,0) 

25116  -  40754 

OVERLAY  (KTIRE,5,0) 

25116  -  27002 

OVERLAY  (KTIRE,6,0) 

25116  -  27034 

Since  overlays  0  to  3  deal  with  input  data  processing  only,  the 
actual  field  length  requirements  for  execution  will  be  estimated  from  the 
blank  cormion  length  calculations  for  overlays  4  to  6. 

Blank  common  length  for  overlay  4  =  24  x  NOEL  +  18  x  NNODE 
+  INCR  X  2  +  MAXLAY  x  10 

Blank  common  length  for  overlay  5  =  3  x  NOEL  +  12  x  NNODE 
+  NEL  +  MAXD 

Blank  common  length  for  overlay  6  =  NEL  +  2  X  MAXD 

where 

NOEL  =  number  of  elements, 

NNODE  =  number  of  nodes, 

INCR  =  number  of  increments, 

MAXLAY  =  maximum  number  of  layers  per  element, 

NRES  =  number  of  restrained  generalized  displacements. 
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MAXD  -  5  X  NNODE  -  NRES,  and 
NEL  =  (MAXD  X  (MAXD  +  l))/2. 

Based  on  the  above  schedule,  the  following  table  provides  sufficient 
data  to  establish  the  field  length  requirements  for  execution  for  a  rota¬ 
tional  ly  symmetric  problem.  Note  that  the  current  version  of  the  code  em¬ 
ploys  in-core  processing  in  the  fifth  and  sixth  overlays.  Approximately 
ll.OOOg  should  be  added  to  these  figures  on  the  Scope  3.4.3  system. 

At  this  time,  only  limited  data  are  available  to  accurately  esti¬ 
mate  execution  times.  Approximately  250  CP  seconds  were  used  on  the  CDC- 
6600  machine  under  Scope  3.4.3  to  perform  8  increments  for  the  following 
problem  size: 

Number  of  elements  *  80, 

Number  of  nodes  =  62,  and 

Number  of  layers  =  6. 

Thus,  as  a  preliminary  estimate,  one  may  consider  4  CP  seconds  per  iteration 
for  10  elements.  The  10  time  is  roughly  twice  the  CP  time  requirement. 
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SECTION  VII 


OPTIONAL  OUTPUT  AND  PLOT  GENERATION 


At 

on  TAPE 2 

1. 

2. 

3. 

4. 

5. 

6. 


each  incremental  step,  the  following  data  are  written 
of  the  TIRE  code: 

SOLUTION  VECTOR  DUE  TO  INFLATION  AND/OR  ROTATION 

STRESS  RESULTANTS 

REFERENCE  SURFACE  STRAINS 

STRESS  DISTRIBUTION 

PRESSURE  INCREMENTS 

CROWN  DISPLACEMENT 


The  above  data  are  then  used  in  an  independent  computer 
code  called  WRITE,  which  may  be  used  to  select  the  desired  out¬ 
put  and  plotting  routine  according  to  the  following  schedule. 

The  first  input  card  contains  the  overall  characteristics 
of  the  problem; 

NOEL  NNODE  NOLAY  INCR  NCROWN 


where  NCROWN  is  the  crown  node  number. 

The  next  card  contains  the  print  and  plot  option  keys; 
KREF  KSTRES  KDISP  KCROWN  KPRES  KCRPLT  KSELCT 

where 


KREF 

KSTRES 

KDISP 

KCROWN 


1;  Provides  resultants  and  reference  strains  for 

each  element. 

1:  Provides  layer  stresses  for  each  element. 

1;  Provides  nodal  displacements. 

1;  Provides  crown  displacement. 
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KPRES  = 

1: 

Provides  pressure  increments. 

KCRPLT  = 

1: 

Crown  displacement  is  plotted  versus  in¬ 
flation  pressure. 

KSELCT  = 

0: 

Stress  resultants.  Strains  and  Layer 

Stresses  are  provided  at  all  incremental 

steps . 

KSELCT  = 

N: 

Stress  resultants.  Strains  and  Layer  Stresses 
are  provided  only  at  the  N  increment. 

Any  combinations  of  the  above  keys  may  be  activated. 

Input  zero  if  a  particular  option  is  not  desired,  with  the  excep¬ 
tion  of  the  KSELCT  key. 
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